Excited states of the helium-antihydrogen system. 
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Potential energy curves for excited leptonic states of the helium-antihydrogen system are cal- 
culated within Ritz' variational approach. An explicitly correlated ansatz for the leptonic wave 
function is employed describing accurately the motion of the leptons (two electrons and positron) in 
the field of the helium nucleus and of the antiproton with arbitrary orbital angular momentum pro- 
jection A onto the internuclear axis. Results for A=0, 1 and 30 are presented. For quasibound states 
with large values of A and rotational quantum numbers J > A no annihilation and rearrangement 
decay channels occur, i. e. they are metastable. 

PACS numbers: 31.30 Jv, 12.20 Ds, 31.15.-p 



Two groups at CERN [l], Q reported recently on the 
production of antihydrogen (H) atoms, which under the 
present experimental conditions were produced in Ryd- 
berg states. Future experiments with fi are aiming for 
investigations of their spectroscopical properties inside 
of a trap. The experimental progress achieved has stim- 
ulated theoretical investigations of atom-antiatom sys- 
tems. Corresponding calculations are usually based on 
variational methods utilizing the Born-Oppenheimer ap- 
proximation. Several works devoted to the H-H system 
were considering both scattering phenomena (see, e.g. 
Refs. [1, Hj]) and the formation of quasibound states 
d, The potential energy curve of the HH system 
was found to be a monotonic function of the internu- 
clear distance R. There exists a so-called critical inter- 
nuclear distance R c , where the electron-positron (e~e + ) 
pair becomes unbound to the nuclei (proton and antipro- 
ton, p + p~). Close to R c the Born-Oppenheimer approx- 
imation breaks down. The optical potential method de- 
veloped in [i[ is aiming for an adequate description of the 
positronium (Ps atom) ejection within the framework of 
the adiabatic picture. The interaction between the He 
and H atom both being in their ground states seems also 
well understood M, UQ. A small potential wall has 
been discovered [9j, which allows for a few quasibound 
states with lifetimes of about 10" 12 - 10" 10 sec 7]. The 
Born-Oppenheimer approach also applies for describing 
the ground state of the HeH system. For investigations 
of antihydrogen-atom collisions non-adiabatic methods 
have been developed as well. Cross sections for scattering 
of H on the hydrogen helium [ll[ and alkali-metal 
atoms 12] have been calculated employing the atomic 
orbital expansion technique. 

Since the H atoms are produced in Rydberg states, 
the atom-H system will be most likely formed in states 
with high values of orbital angular momentum projection 
A onto the internuclear axis. From a pragmatic point 
of view it is important to develop an adequate descrip- 
tion for such states. This was already achieved for the 
HH system in [pi [l4| by utilizing explicitly correlated 



wave functions, which describe accurately the motion of 
leptons with arbitrary orbital angular momentum pro- 
jection. In the present Letter we proceed in a similar 
way and extend our approach to describe the interaction 
between a He atom in its ground (singlet) state and an 
H* atom in an excited state. Atomic units will be used 
throughout. 

Application of the adiabatic approach to the HeH sys- 
tem leads to the Schrodinger equation 



ff lep 4^ p (x,R) = v| op (i?)1^ p (x,R) 



(1) 



for the wave function \l/^ p describing the motion of the 
leptons in the field of the He nucleus (a-particle) and of 
the antiproton Here x = (r l7 r 2 ,r 3 ) denotes the position 
vectors of the leptons and R = R Q — R p defines the 
internuclear distance R = |R| between the nuclei located 
at R Q and R p _ with respect to the center-of-mass frame. 
The leptonic Hamiltonian reads 
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where a — =pl refers to the charge of e~ and e + , re- 
spectively. The HeH interaction energy is expressed in 
terms of the leptonic potential E A (R) = V^ P (R) — 2/R. 
Axial symmetry implies the eigenvalue equation 

L R ^ p (x,R) = A^ p (x,R) (3) 

for the component L R of the leptonic orbital angular mo- 
mentum operator along the internuclear axis with eigen- 
value (orbital angular momentum projection) A. If the 
internuclear distance tends to infinity the He and H 
atoms no longer interact. Thus, in the limit R — > oo 
the leptonic wave function equals the product 



^ ep (HeH) = * gs (He)* es (H) 



(4) 
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between the wave functions \& (He) and \I/ CS (H) describ- 
ing the He atom in the ground state (gs) and the H* 
atom in an excited state (es), respectively. In the case of 
large internuclear distances the leptonic orbital angular 
momentum projection A of the Hefi* system is mainly 
carried by the c + . 

For solving the Schrodinger equation {1} an explicitly 
correlated ansatz for the leptonic wave function is taken 
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where we defined the vectors Vfe = J2 i=1 u^r.^ The coef- 
ficients C k are linear and &y, a*, and are nonlin- 
ear variational parameters, while the operator P ensures 
proper symmetry. In our calculation a set of JV = 300 
basis functions {t/'aIjcLi was adopted. Each function tj}\ 
appears as product of Explicitly Correlated Gaussians 
(ECGs) together with an angular part r/ A . The ECGs 
alow for an adequate description of the lepton in the 
field of the nuclei. The angular part r/ A involving a usual 
spherical harmonic Y AA ensures the proper angular sym- 
metry of each basis function tp A according to equation 

For any fixed internuclear distance each nonlinear pa- 
rameter was optimized employing the golden section 
method. Though there are more refined optimization 
methods, e. g. (l5j . we prefer to use this very simple 
but reliable approach. For few-particle systems such as 
HH or HeH quasimolecules it does not look much more 
cumbersome than other more sophisticated approaches. 
The set of linear variational parameters {C k } k=1 was ob- 
tained by solving the generalized eigenvalue problem. 

Results for the HeH interaction energies E A for A=0, 
1 and 30 are presented in Table 1. The behavior of the 
potential curve for the ground state of the HeH system 
(A=0) was first examined in [9]. Our results for A=0 
are in a good agreement with previous investigations @ 
and achieved with half the number of basis functions. 
Furthermore, we shall focus on the excited levels of the 
HeH system. 

The leptonic potentials V^ ep for HeH with A=0, 1 and 
30 are depicted in Fig. Q] together with the potential Vp s 
for a compound Hc + p~ + unbound Ps system in its low- 
est state. We assume, that for a state with a certain 
value of A the Ps ejection occurs, if the potential en- 
ergy curve V A P coincides with the curve V Ps at distances 
R < B, c (A) . In principle, the crossing of the two curves 
does not immediately imply the ejection of a Ps atom. 
With respect to the nuclear center-of-mass system the 



Table 1: The interaction energy E A (R) as a func- 
tion of the internuclear distance R is calculated for 
leptonic orbital angular momentum projections 
A=0, 1, 30 (in atomic units). 
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FIG. 1: The leptonic potentials for the ground state (A = 0) 
and for excited states (A = 1, 30), Vg Cp (R) (dash-dotted line), 
Vl ep (R) (solid line) and V^ P (R) (dashed line) are plotted 
versus the internuclear distance R (in atomic units). The 
potential for the lowest state of the He + p~ -(-unbound Ps 
system V Pa (R) (dotted line) is also plotted. The potential 
V 1 lep (i?) approaches V Ps (R) near the critical internuclear dis- 
tance R c — 0.7a . 
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two parts: the binding energy £ T ,1 s Ild of the ground state 
of the Ps atom and the rotational energy Pp°* of the Ps 
atom with respect to the intcrnuclcar axis. For values 
A > but small, the energy E P °^ is negligible compared 
to |-Ep s nd |. Conversely, for A > 1 the contribution 
is essential as it was shown in [3, 0] for the HH* sys- 
tem. The rotational part of the energy Ep°^ is not taken 
into account in the potential for the Hc + p~+unbound Ps 
system, since it is inversely proportional to the square of 
the distance between the Ps atom and the He + p~ com- 
pound. Thus, the crossing of the potential curves V Ps and 
V^ cp does not mean that the energies of the two systems 
(HeH* and He + p _ +Ps) become equal. For the states 
with large values of A this takes place only, when the in- 
ternuclear distance R becomes sufficiently small and the 
Ps atom appears to be far enough from the He + p~ com- 
pound, so that the rotational energy becomes negligible. 
Consequently, the curves for the HeH* and the combined 
He + p~+Ps systems coincide for all values of R < R c . 

According to Fig. [1] the properties of the states of the 
HeH system with the A values under consideration (A=0, 
1, 30) differ strongly. The HeH potential energy curve 
Vq CP does not cross the one for the He + p~+unbound 
Ps system. Therefore, the light particles appear to be 
bound in the HeH system over the entire range of inter- 
nuclear distances R. For R — > the function V^ lep ap- 
proaches the binding energy of the positronium hydride 



E RPs = -0.7891967 [16J. This reveals that at small in- 
ternuclear distances the HeH quasimolecule transforms 
into the He + p~ system plus the Ps atom weakly at- 
tached to it. The HeH* states with nonzero but small 
A (see e.g. A=l) exhibit properties similar to those of 
the Hfi* system [H, There exists a critical internu- 
clear distance i? c (A = 1) = 0.7a at which the wave func- 
tion \l/^ p transforms from the wave function of a bound 
HeH* system into that of an unbound Ps atom in the 
field of the He + p~. The latter wave function then con- 
tains a plane wave factor describing the center-of-mass 
motion of the Ps atom. Thus, the adiabatic correction 
to the leptonic potential diverges near R c indicating the 
breakdown of the Born-Oppenheimer approximation in 
the vicinity of the critical distance. Despite the fact that 
the ECGs cannot properly reproduce a plane wave, we 
keep the ansatz Eqs. (|5I6I7[) even for internuclear dis- 
tances R < R c . However, the basis set (JBJ) provides the 
correct value for the energy of the lowest continuum state 
of the He + p _ -{-unbound Ps system (i. e. with zero rela- 
tive velocity of Ps atom with respect to He + p~ ) as well as 
for the matrix element involving a spatial delta-function 
(see below). Properties of the states with large values 
of A have already been elucidated for the HH system in 
[HI, [HI- As implied by Eq. Q, for large internuclear dis- 
tances the HeH* system can be envisaged as a He atom 
plus a p~ and e + weakly attached to them. As the inter- 
nuclear distance R decreases the system changes slightly 
in the following sense: The orbital angular momentum 
can still be mainly attributed to the e + . Actually, the 




FIG. 2: The coalescence probability distributions Po(R) 
(dash-dotted line), Pi{R) (solid line) and P 30 (R) (dashed line) 
are plotted as a function of the internuclear distance R (in 
atomic units). 



potential curves for the HeH* and Hep systems behave 
similarly [l7|. At R = the potential V^q P approaches 
the value of the binding energy of an H~ ion plus a small 
contribution from a weakly bound e + (see Fig. [I}. Hfi* 
and HeH* systems in states with high A, respectively, 
differ mainly in their behavior at small internuclear dis- 
tances. When the p~ approaches the p + in Hfi* system, 
the e~ becomes loosely bound and forms a Ps atom to- 
gether with the e + . In case of the HeH* system, the 
a-particle tends to keep both e^s due to its larger elec- 
tric charge and the Ps atom is not ejected. According to 
Fig. [I] the leptonic potentials for the Hefi* with A = 30 
and that of the He + p~+unbound Ps system cross each 
other. However, as explained above, this does not imply 
the ejection of the Ps atom for the HeH* system in the 
state with A = 30. 

The results discussed above can be confirmed by eval- 
uating the coalescence probability distribution P^ 7 Q as 
a function of R, which arises in calculations of the two- 
photon leptonic annihilation rate. According to [3] the 
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general expression for P A is 
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where the spin operator 1 — Sf j ensures that the 2j- 
annihilation can only take place between c~ and c + be- 
ing in the singlet state. The annihilation rate of a bound 
state of a particle and antiparticle is proportional to the 
matrix element of the spatial delta function. The coales- 
cence probability distributions for the HeH system with 
A=0, 1 and 30 are depicted in Fig. [21 If the Ps atom 
is ejected, the function P A 7 should approach the value 
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of the coalescence probability for the ground state of the 
positronium IVW )! 2 = 1 /§ 7r ~ °- 038 - 

In the highly excited HeH system with A=30 the den- 
sity of the very weakly bound e + is thinly distributed 
over the entire range of internuclear distances. Thus, 
the probability for the e + to coalesce with the e~ is 
negligible for all values of R (see Fig. [5]). On the other 
hand, for A=30 the Ps atom (bound or unbound) is never 
formed within this adiabatic picture. Consequently, the 
lcptonic annihilation decay rates and the Ps ejection de- 
cay rates are negligible for the HeH system in the states 
with high values of A. In particular, for the quasibound 
states of HeH* with A=30 and a rotational quantum 
number J > A the angular momentum barrier prevents 
the a-particle and p - to coalesce. Hence, the nuclear an- 
nihilation decay rates vanish for such quasibound states. 
Since there are no annihilation and Ps ejection decay 
channels, the Hcfl* quasimolecule becomes metastable in 
these states. The decay of such states with high values of 
A and J > A can occur only via a radiative cascade into 
a final state, where the annihilation process is probable. 
The lifetime of such metastable states of the HeH* quasi- 
molecule are expected to be of order 10 -6 s as for Hep~ 
atomcules [l9l |20| ; as mentioned above the properties of 
these systems are similar. 

Under the experimental conditions reported in [l|, Q 
the H atoms are produced inside of traps with very high 
magnetic fields (up to 5 T). Atomic levels with differ- 
ent angular quantum numbers I will be fully admixed 
and remain so, when the quasimolecules are formed. If 
the molecular axis is oriented parallel to the magnetic 
field, its presence will not lead to any qualitative differ- 
ence compared to the case of zero field. According to the 
"guiding center atoms" picture of the three-body inter- 
action in a plasma developed in [2l|, [12] , the formation 
of HH and HeH quasimolecules should be most proba- 



ble with this orientation. Moreover, the magnetic field 
influences more strongly the quasimolecular (rotational) 
levels rather than the behavior of the potential curves; 
the latter being the major subject of the present Letter. 

Summarizing, we can state that the method employ- 
ing an explicitly correlated ansatz as developed recently 
for describing Rydberg states of the Hfl system, can be 
successfully applied for calculations of the HeH system. 
Accurate potential energy curves for the quasimolecular 
states with A=0, I and 30 are obtained. The results for 
the ground state of the HeH system (A=0) are in agree- 
ment with the known ones. The potential energy curves 
and coalescence probability distributions obtained can be 
used for evaluating of the spectrum of the HeH quasi- 
molecule, decay rates and cross sections for various pro- 
cesses. The prediction of metastable states in the HeH 
system (e.g. A=30 and J > A) is the most important 
result of the present Letter. This leads to the possibil- 
ity to deposit a big amount of energy at atomic scales 
(~ I GeV per molecule) over a relatively long time pe- 
riod (~ 10~ 6 s). Assuming typical atom velocities (~ I0 7 
cm/s) this implies a possible transfer of this energy over 
macroscopic distances. The problem of coexistence of 
matter and antimatter is of fundamental interest not only 
for laboratory studies but also for cosmology. The enor- 
mous enhancement of the lifetime of Rydberg states of 
the HeH* quasimolecule compared to HH* (about I0 8 ) 
may have most important consequences. 
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